Corotation: its influence on the chemical abundance pattern of 

the Galaxy 
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ABSTRACT 



A simple theory for the chemical enrichment of the Galaxy which takes into account the effects 
of spiral arms on heavy elements output was developed. In the framework of the model with the 
corotation close to the position of the Sun in the Galaxy the observed abundance features are 
explained. 

Subject headings: Galaxy: abundances - Galaxy: spiral arms: Galaxy - corotation 



1. Introduction 

The radial distribution of heavy elements in 
the Galaxy is a signature of the Galactic history 
and structure. Especially important information 
is brought by various features in the abundance 
pattern. 

In early papers the radial variation of abun- 
dance was approximated by a function close to the 
linear one in logarithmic scale (see e.g. Harris, 
1981; Harris & Pilachowski, 1984 and references 
therein). However, in the last decade some fea- 
tures of the chemical distribution were revealed 
that do not support the simple linear model. In- 
deed, Simpson et al (1995) and Vilchez & Esteban 
(1996) derived a plateau- like (or step- like) abun- 
dance structure of the H II regions. The same 
pattern was found by Twarog et al (1997) for open 
clusters, by Smart & Rollestone (1997) for B stars, 
by Caputo et al. (2001) for Cepheids. In recent 
papers Andrievsky et al. (2002, a,b) on the basis of 
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precise spectroscopic studies of Cepheids derived 
a "bimodal" abundance gradient with a plateau 
beyond the solar distance and a sufficiently steep 
gradient inside the solar circle. Finally, Maciel & 
Quireza (2000) found a parabolic metallicity dis- 
tribution for the planetary nebulae, with a mini- 
mum near the solar distance. 

What is the mechanism of such features forma- 
tion? In this paper we will inspect the idea that 
the bimodal abundance gradient can be a conse- 
quence of the corotation resonance. The underlying 
idea is that at corotation the star formation rate 
presents a minimum, and consequently, the metal- 
licity accumulated after a few billion years at this 
radius can also have minimum. Bearing in mind 
that since Marochnik et al (1972) and Creze & 
Mennessier (1973) a lot of studies showed that the 
corotation is situated very close to the Sun (Nelson 
& Matsuda, 1976; Mishurov et al, 1979; Mishurov 
et al, 1997, Amaral & Lepine, 1997; von Linden 
et al., 1998; Leicht & Vasisht, 1998; Mishurov & 
Zenina, 1999a, b; Lepine et al, 2001; Fernandez et 
al., 2001 etc) it would be interesting to analyze 
the influence of the corotation on the abundance 
pattern. 

On the other hand, according to Weinberg 
(1994), Englmaier & Gerhard (1999), Dehnen 
(2000) and others the corotation is situated very 
close to the galactic center at a distance ~ 4 kpc. 
So, investigating the chemical abundance struc- 
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ture of the Galactic disk, we can hope to make a 
choice between the above models. 



a form close to that of LF: 



2. Basic idea and equations 

Standard equations governing the chemical evo- 
lution of the galactic disk take into account both 
the creation of heavy elements and mass transfer 
between gaseous and stellar components (Tinsley, 
1980). They are essentially the equations of conti- 
nuity and convey the mass conservation law for a 
mixture: gas + stars + heavy elements (e.g. Lacey 
& Fall, 1989, hereafter LF; Portinari & Chiosi, 
1999, hereafter PC). By means of these equations 
researchers usually try to simultaneously explain 
both the radial behavior of metal abundance and 
gaseous and stellar distributions. However the cor- 
responding equations must be completed by the 
dynamical ones for instance by the Euler gasdy- 
namic equation, the Boltzman equation for the 
stellar component, the Poisson equation for grav- 
itation field etc. Only such an approach would 
enable us to derive all radial profiles in a self- 
consistent manner. Of course, this way is very 
complicated and most authors prefer to exclude 
the dynamical part of the problem (a rare excep- 
tion is the paper of Martinet & Friedly, 1997) and 
consider the task in the framework of the stan- 
dard scheme, specifying some model laws for the 
star formation rate (SFR), the infall rate of matter 
onto the Galactic disk, and the radial flows within 
the disk, so as to explain the observed radial dis- 
tribution for chemical abundance, star density and 
interstellar gas density. 

To make more prominent the effects of the spi- 
ral arms and of the corotation on the abundance 
pattern which are not masked by other factors, in 
the present paper we shall only concentrate our at- 
tention on the chemical enrichment. For this task 
we adopt the stationary model of Haywood et al. 
(1997) for Galactic evolution. According to these 
authors, gas density and infall rate are constant 
in time. This enables us to consider the chemi- 
cal part of the problem separately for various laws 
of infall rate, radial flow and star formation rate 
(SFR) 4 . The corresponding equation we write in 



dZ 

dt 



4 Of course, we abstain from explaining the stellar density 
distribution. In any case it is impossible to adequately 
derive the complicated structure of the stellar disk in the 
framework of models which do not take into account dy- 
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where \i g is the surface density of interstellar gas, 
Z is the abundance of heavy elements, ip is the 
star formation rate, Pz is the mass fraction of a 
star ejected as newly synthesized heavy elements 
(Tinsley, 1980), u is the macroscopic radial veloc- 
ity of the gas, / is the infall rate of matter onto the 
Galactic disk per unit area, Zf is the abundance 
of infalling gas, t is time, R is galactocentric dis- 
tance. In addition to the corresponding equation 
of LF we include the effect of diffusion (in the last 
term in the right hand side of Eq (1), D is the 
coefficient of heavy elements diffusion). The mat- 
ter is that, due to stellar winds, supcrnovae ex- 
plosions, Galactic fountains, etc., the interstellar 
gas is highly turbulent. This leads to gas mixing, 
which tends to smooth out the abundance gradi- 
ents. For a phenomcnological description of this 
process we introduced the diffusion term into Eq. 
(!)• 

From previous studies it is known that the in- 
fluence of the SFR radial law happens to be very 
important on the Galactic abundance structure. 
Numerous approximation formulas were proposed 
for it. As a rule the SFR ijj is assumed to be pro- 
portional to some power of gas surface density /j, g . 
However this cannot be recognized as completely 
satisfactory. Indeed, it is well known that super- 
novae of all types - sources of heavy elements - are 
strongly concentrated in spiral arms (Bartunov et 
al, 1994). Hence synthesis of heavy elements in 
a given volume of matter takes only place when 
the volume occurs to be inside an arm. To take 
into account the effects of spiral arms and corota- 
tion on chemical enrichment we represent the SFR 
in a manner similar to Oort (1974), Wyse & Silk 
(1989) and PC: 



V = /3^|0(fl)-n P |erv) 



(2) 



where Q(i?) is the angular rotation velocity of the 
Galactic disk, ftp is the rotation velocity of the 



namical processes like Galactic contraction as a whole, in- 
fluence of the bar and of the spiral structure etc. 
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Galactic density waves responsible for spiral arms 
(recall that the spiral pattern rotates as a solid 
body, i.e. Qp = const, whereas the Galactic disk 
rotates differentially, i.e. f2 is a function of R, the 
distance Rc where both the velocities coincide, 
i.e. Vl(Rc) = $lp, is called the corotation radius), 
(3 and k are some constants (we select the product 
Pzf3 so that Z(Rq) = Zq, 5 and consider two cases 
k = 1 and k = 1.4 Kennicutt, 1998), 6 is a cutoff 
factor (see below). 

The factor — ttp\ in Eq (2) means a frequency 
of crossing spiral arms by some element of matter 
of Galactic disk. It appears due to the fact that 
the enrichment by heavy elements takes only place 
when the element of Galactic matter occurs to be 
within a spiral arms. An equivalent interpretation 
is that the SFR depends on the rate at which the 
gas penetrates the spiral arms. 

Wyse & Silk (1989) and PC considered the case 
when the corotation resonance is situated at the 
very end of the Galaxy (model of Lin ct al, 1969). 
So, they simply neglected by the value fip in Eq 
(2). 

In our paper we do not neglect fl p and consider 
several cases: with the corotation close to the po- 
sition of the Sun (Rc ~ Rq = 8.5 kpc); corotation 
near the galactic center (Rc = 4 kpc) and far coro- 
tation (Rq = 14 kpc). 

The cutoff factor is required because the self- 
sustained Galactic density waves are known to ex- 
ist only in the so-called wave zone between in- 
ner Lindblad resonance i?j„ where v(Ri n ) = — 1 
and the outer Lindblad resonance R out where 
v(Rout) = +1 (Lin ct al., 1969). So, we de- 
fine 6 = 1 for \u\ < 1 and 6 = for \v\ > 
1. Here, v = m(Clp — £Y)/k is the dimension- 
less wave frequency, m is the number of arms, 
k = 20(1 + Rdn/dR/2Q)^ 2 is the epicyclic fre- 
quency. 

The expression (2) for SFR means that we only 
consider an ordered star formation process trig- 
gered by spiral arms. Besides this, in the Galaxy, 
some sporadic (or spontaneous) star formation 
may take place. Hence, in general ip should include 
the corresponding term. But we restrict ourselves 



5 According to various authors (see e.g. PC and references 
therein) the solar metallicity is greater than that of its sur- 
roundings. However Andricvsky et al. (2002a) showed that 
the Sun has a metallicity which is typical of its vicinity. 



to the representation (2) for the SFR; in other 
words, we consider that the spontaneous SFR is 
negligible compared to that induced by the spiral 
arms. 

From the above equation, for a given set of in- 
put functions, we can compute the time history 
of chemical evolution of the Galactic disk. How- 
ever, in this paper we only consider the simplest 
asymptotic solution. Indeed, due to diffusion term 
in Eq (1), we expect that dZ/dt — > for t — > oo. 
A typical diffusion time for the achievement of the 
asymptotic state is tp, <~ L 2 /D where L is a typ- 
ical length, say L ~ 1 kpc. Let us estimate tp> 
and show that we may use the asymptotic regime. 
For our phcnomcnological theory we consider in- 
terstellar gas clouds as particles. So, we use the 
gas kinetic theory for the estimation of the diffu- 
sion coefficient: 



a ; 



H g D 



m c 



HV T 



where mc is the mass of a cloud, a is its typical ra- 
dius, H is the thickness of gas layer (H = 0.13 kpc) 
and Vp is the thermal velocity of a cloud (Vp = 
6.6 km s^ 1 , Stark, 1984). Following Falgarone & 
Puget (1986; see also Elmegreen, 1987) wc adopt 
mc/a 2 = 100Mq/pc 2 . Notice that for constant 
H and Vp the value a is constant as well. For a 
typical /j, g ~ 10Mq/pc 2 tp> <~ 1.6 Gyr. So, in the 
major part of the Galaxy we may use the above 
asymptotic approximation although in a more re- 
fined theory departures from it should be investi- 
gated. 

As boundary conditions we adopt: at the "end" 
of the galactic disk, say at R = Rc — 25 kpc 
Z = Z f (Z f = 0.17 Z Q , Z q = 0.02 Tinsley, 1980) 6 
and at the galactic center the abundance must be 
finite, therefore dZ/dR = for R = (in what 
follows we shall not be interested in the abun- 
dance structure of the central part of the Galaxy 
since chemical composition there is masked by 
very complicated processes due to the bar, the 
bulge etc.). 

Let us now discuss the functions entering Eq 
(1). For the infall rate we adopt: 



/ = foexp(- 



R — Rq 



0- 



3 Results depend weakly on the exact value Rq in the range 
20 - 30 kpc 
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where /© = 3.5 M Q pc 2 Gyr 1 (Haywood et al, 
1997). 

For the surface gas density we use a simple rep- 
resentation 

, R d, 
fig = n Q exp{-- - -) 

which imitates the hole in the central part of the 
disk with maximum n g at R = d and is close 
to observed radial distribution of interstellar gas 
(Dame, 1993) for d = 4 kpc, the dependence be- 
ing normalized as to give /x 9 (i? Q ) = 10 M Q pc~ 2 
(Haywood et al, 1997). 

In all computations the rotation curve of 
Allen & Santillan (1994) with R Q = 8.5 kpc and 
V ro t{RQ) — 220 km s^ 1 was used. 

3. Results and Discussion 

First of all we consider a model with corota- 
tion close to the Sun (Rc — Rq) and radial inflow 
u = 0. Our experiments with A in the range 4-5 
kpc (close to usually adopted) give similar results. 
In Fig 1 by solid lines are shown the abundance 
patterns for k — 1.4 (thick line) and k = 1 (thin 
line) for A = 5 kpc. To feel the effects of coro- 
tation, diffusion and cuttoff 6-tcrm on chemical 
enrichment we give the stationary solution of Eq 
(1) neglecting the diffusion (k = 1.4 - dotted line; 
k = 1 - dashed line). Note abrupt decrease of 
abundance in the last two cases beyond the Lind- 
blad resonances due to cuttoff 6-term and the gap 
in the vicinity of the corotation radius. As it was 
expected those sharp features are smoothed out 
due to diffusion. 

In Fig 2 thick line is the same as in Fig 1 but 
the dashed line represent the solution of Eq (1) for 
the corotation slightly closer to the galactic cen- 
ter, namely Rc = 8kpc and k = 1.4. Both the 
curves are superimposed on the observed abun- 
dance derived from Cepheids by Andrievsky et al. 
(2002 b). From this Figure it can be seen that 
the model with the corotation close to the Sun 
describes rather well the bimodal abundance gra- 
dient derived by Andrievsky et al (2002 b), with 
a plateau beyond the solar circle and a steep gra- 
dient inside of it. 

In Fig 3 is shown the influence of radial gas 
flow: i) with constant velocity u: u = +0.5 kms^ 1 
and u = — 0.5 kms^ 1 and ii) with outflow from 



the corotation with the same amplitude 0.5 kms 1 
(Lepine et al, 2001). Thick line is the same as in 
Fig 1. 

From the above Fig 1-3 one can see that it is 
possible to explain the abundance features like bi- 
modal abundance structure of Andrievsky et al 
(2002) by means of combined influence of spiral 
arms, closeness of corotation to the solar galacto- 
centric distance and radial gas flow on chemical 
enrichment. 

In Figs 4,5 by solid lines the results are given 
for models with the corotation near the galactic 
center (Rc — 4 kpc, Fig 4) and far corotation 
(Rc = 14 kpc, Fig 5). In both pictures dotted lines 
represent stationary solutions obtained by neglect- 
ing the diffusion term. From these picture it is ob- 
vious that these models for Galactic spiral waves 
are not consistent with the observed abundance 
pattern. 

We are grateful to S.M. Andrievsky for very use- 
ful discussion of the problem. 
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Figure captions 

Fig. 1 - The theoretical radial abundance distri- 
bution for the corotation radius close to the Sun 
position, i.e. Rc — Rq — 8.5 kpc. Solid lines: 
thick - k = 1.4, thin - k = 1. Dotted and dashed 
lines describe the stationary solution of Eq (1) ne- 
glecting the diffusion term for k = 1.4 and k = 1 
correspondingly. Radial inflow u = 0. 

Fig. 2 - The theoretical abundance gradient 
(solid and dashed lines) superimposed on the ob- 
servational data derived by Andrievsky et al (2002 
b) for Ccphcids (crosses). The thick line is the 
same as in Fig 1, the dashed one corresponds to 
the corotation radius slightly closer to the Galac- 
tic center then in the previous case, namely Rc = 
8 kpc. Radial inflow u = 0. 

Fig. 3 - The theoretical abundance gradient 
taking into account the radial inflow: i) u = 
+0.5 ferns -1 - dashed line; ii) u = — 0.5 kms^ 1 

- dotted line; iii) outflow from the corotation 
(Lepine et al, 2001) with the amplitude 0.5 kms^ 1 

- dashed - dotted line, the corotation radius being 
situated at the Sun galactocentric distance. The 
thick solid line is the same as in Fig. 1. 

Fig. 4 - The theoretical abundance gradient 
for the model with the corotation near the galac- 
tic center, Rc — 4 kpc and k = 1.4 - solid line. 
Dashed line represents the stationary solution ne- 
glecting the diffusion. 

Fig. 5 - The same as Fig. 4 but for the far 
corotation model, Rc = 14 kpc. 
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